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Abstract 

The present work studies a kind of Maximum Concurrent Flow Problem, 
called as Extended Maximum Concurrent Flow Problem with Saturated Ca- 
pacity. Our major contributions are as follows: (A) Propose the definition 
of Extensive Maximum Concurrent Flow Problem with Saturated Capacity 
and prove its solutions exist. (B) Design a approximation algorithm to solve 
the problem. (C) Propose and prove the complexity and the approximation 
measures of the algorithm we design. 

Keywords: network; concurrent flow; approximation algorithm; 
approximation measure. 
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1. Introduction 

The Maximum Concurrent Flow Problem (MCFP) is a main kind of mul- 
ticommodity flow problems, which was introduced by Matula in 1985, see 
Shahrokhi and Matula [9], and have been severely studied for more than two 
decades, see, e.g., [2,4,5,6,9,10]. Motivated by the pubhshed literatures of 



Preprint submitted to Elsevier 



June 28, 2011 



this research area, we discuss a kind of specific Maximum Concurrent Flow 
Problem in the present work, termed as Extended Maximum Concurrent 
Flow Problem with Saturated Capacity (EMCFPSC). 

The rest of this article is organized as follows. Some preliminaries are 
presented in Section 2. Section 3 formulates the problem EMCFPSC and 
prove the existence of its solutions. Section 4 is specially devoted to designing 
an algorithm to approximately solve the problem. Section 5 investigates the 
complexity and precision of the algorithm. Finally, the paper is concluded 
with Section 6. 

2. Preliminaries 

This section provides some preliminaries for our sequel research. 

A graph G — {V,E) is called as a hybrid graph if its edge can be either 
directed edge or undirected edge, of which the directed graph and undirected 
graph are specific cases. Here V and E represent all the nodes (or vertices) 
and edges (or arcs) of G respectively. 

Given graph G= {y,E) with weight (or capacity) c and H = {[si,ti] : 
Si,ti & V;i — 1,2, - ■ ■ ,k}, where c : E — > i?+ and R+ — [0, oo), and 
Si and ti express the source and the terminal of commodity i respectively, 
we call the triad {G,c,H) as a multicommodity network. Let [si,ti] e H; 
Si,Vi,V2, ■ ' ' I'^hti^^ some nodes of G and be different from each other except 
for Si = ti] Cj be the edge of G with endpoint Vj resp. Vj+i for < j < I, 
where Vq — Sj and vi+i — ti. (When is directed, the edges Vj and Vj+i must 
be the original endpoint and the terminal endpoint, respectively.) Then we 
call P = [{si, eo,vi), {vi,ei,V2), ■■ ■ , {vi,ei,ti)] as an Sj - path of {G,c,H). 
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(To simplify in notation afterwards, we denote {vj, ej, Vj+i) as {vj, Vj^i) when 
Cj needn't be indicated; and denote {vj,ej,Vj+i) as ej when vj and vj+i 
needn't be indicated.) Let Vi be a set of Si — ti— paths such that c(e) > 
for aU e e E{P) and P e Vi, i = 1,2, - ■ ■ ,k, where E{P) denotes the set 

k 

of all the edges of P. We call V = |J T^j as a (positive) path system on 

i=l 

{G, c, H), which is denoted by V\(g,c,h): or V for conciseness. 

Definition 1. Let V he a path system. If mapping y : V — > R+ satisfies: 

y{P)<c{e) yeeE{V), 

eeE{P),Per 

where E{V) = {e e E{P) : P eV}, i.e. the set of all the edges of V. Then 

we call y as a multicommodity flow on V defined by the function on the path 

system, and as a flow for simplicity; and call V{y) — ^ y{P) as flow value 

Per 

of y, and Vi{y) = ^ j/(P) as branch flow value of commodity i with y. The 

P&Vi 

set of all the flows on V is denoted by F[V], and by F in brief. 
Remark 1. As is known to us, there are two kinds of definitions about 
the flow, which are the defined by the function on the edge set and by the 
function on the path system. For the purpose of being simple and clear, we 
only make the study as far as the flow deflned by the function on the path 
system in the present work. 

Generally, the problem to find a multicommodity fiow to satisfy certain 
conditions is called the multicommodity flow problem (MFP). A few of usual 
multicommodity flow problems are as follows. The problem to flnd a flow 
y so that V{y) is maximized, i.e. V{y) — OPT[P] — V{— max{V{y) : y e 
F[P]}), is called the maximum multicommodity flow problem on V (MMFP). 
A flow y satisfying MMFP is called a solution of MMFP. The problem to find 
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a flow y so that Vi{y) < hi for given bi, i — 1,2, - •• ,k, and V{y) is maximized, 
i.e. V{y) = V{= max{V{y) : Vi{y) < bi,i ^ 1,2, k;y e F}), is called 
the maximum flow problem with bounds(MMFP-B). The problem to find a 
flow y so that V{y) is maximized under the condition Vi{y) = Xbi for given 
hi, i — l,2,---,k, and for A e [0,1] is called the Maximum Concurrent 
Flow Problem (MCFP) (see, e.g., Shahrokhi and Matula [10], or Garg and 
Konemann [6]). 

3. Model formulation 

Now we begin to introduce the problem we are attacking in the present 
work. 

Definition 2. Suppose P is a path system on the multicommodity network 
{G, c, H), and b = {bi,b2, ■ ■ ■ ,bk)- We call the problem to flnd a flow y' so 
that min [r:Vi(?/')] = max{ min [r:V^i(y)] : y e F, Vi{y) < bi} as Extended 

l<i<fc » l<i<fc » 

Maximum Concurrent Flow Problem (EMCFP). We call the problem to find 
a fiow y' so that min [^Vi(|/')] = max{ min [^Vi(i/)] : y E F, Vi{y) < b^}, and 

l<j<fc ' l<i<k * 

V{y') = max{y(y) : y is the solution of EMCFP } as Extended Maximum 
Concurrent Flow Problem with Saturated Capacity (EMCFPSC). 
Theorem 1. The solution of the problem EMCFPSC exists. 
Proof. It is obvious that the solution of MCFP is also the solution of EMCFP. 
Hence, B = {y E F : y is the solution of FMCFP} ^ from the well known 
fact the solution of the problem MCFP exists. When the member of B is 
finite, the conclusion of Theorem 1 is trivial. Otherwise, sup{V{y) : y e B} 
exists and is finite. Let it be b. Then there is at least a sequence {yj} in 

k 

B such that hm [^ V{y) ~ b, and hm yj{P) exists for all P eV. Now let 
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y{P) = lim yj{P) for all P eV. Then y =e B, and V{y) = 6. That is, y is 

i-*-oo 

the solution of EMCFPSC. □ 

Example During the winter between the year 2009 and 2010, in major re- 
gions of the world, the amount of energy consumed for heating is far more 
than the ordinary winter due to the abnormal lower air temperatures. With 
this circumstances, the electrical networks emerge to congestion for a lot of 
districts in the world. Assume the purpose of supplying power is firstly to 
maximize the minimum of the met rates of cities. Then the problem to opti- 
mize the power supplement can be largely regarded as the mould EMCFPSC. 
One can believe the mean of apphcations of the mould EMCFPSC from the 
instance. 

4. Algorithm 

In this section, we are devoted to designing an approximation algorithm 
for the problem EMCFPSC. 

Algorithm (Approximation algorithm for EMCFPSC): 

Input Multicommodity network {G, c, H), path system V, vector b and error 

parameter < < 1. 

Output Flow y* on V, which has the characters specified by the Theorem 
in Section 5. 

1. Put e = min{^, |}; 6° = 6^, z = 1, 2, • • • , fc; / = 0, /i = 0. 

E bi 

i=l 

2. Set I — l + l. Ulr] > 1, implement the next Step. Otherwise, put I — I* 
and go to the Step 4. 

3. Set bi = Itjb^, b = {bi,b2, - ■ ■ , bn). For V, b and e, obtain an approxi- 
mation solution yi of MMFP-B with by Algorithm 2 of Cheng [3]. If 
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^{yi) < TT^(X] ^«)) put I — I* and implement the next step. Otherwise, 

i=l 

return to (2). 

4. Construct the auxihary network and demanding vector b as follows. Set 
V = V\J{t'„t[ : i = l,2,.--,k}{t'„t[eV),E' = E\J{{U,t',),{U,t'^ : 
i = l,2,---,k},H= {[s„ t'l [s„ Q:t = l,2,---,k},G' = (V, E'); V'i = 
{P+{U, t'^:Pe P,}, n = {P+ g : P G P,}, i = 1, 2, ■ ■ ■ , k- V = 
( U V'^ U( U ^'o); = - 1)^6°, = 6, b' = 6^, • • • , and 

i=l 1=1 



c'(e) 



c(e) e e ^(G) 



k 

5. Put h := h + l,b = hr). If 6 > ^ 6°, put h* = h and implement the 
next step. Otherwise for V' , b' and e, obtain an approximation solution 

k 

y2 by Algorithm 2 of Cheng [3]. If V{y2) < i^[(E + P^t /i* = /i 

i=l 

and implement the next step. Otherwise, put y' — y2, and then return 
to 5. 

6. Put y*{P) = y'{P + {U, t'^) + y'{P + {U, to)), VP e V^,^ = 1,2,- ■ ■ ,k. 
Stop. 

5. Algorithm analysis 

This section is specially devoted to discussing the correctness, approxi- 
mate precision and complexity of the algorithms we nave presented above. 
Theorem The complexity of the Algorithm is 0(^A;m^ logn), where k — 
\H\,n ^ \V{G)\,m = \E{G)\. Let y be a solution of EMCFPSC and y* 
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be the output of the Algorithm. Then y* is a flow on Vi{y*) < bi,i — 

k 

l,2,...,k,V{y)<[l*{Zh) + h*]rj,8^nd 

i=l 

k k 

[(r - k) + {h* - l)]rj -2r]< V{y*) < [{1* - k) + h*]rj. 

i=l i=l 

n* u 2r^ ^ . V,{y*) . V,{y) 

[I — 1)7] — < mm — ; < mm — ; — < I rj. 

min bi bi bi 

Proof. It is obvious that the complexity of the Algorithm depends on the 
complexity of the subroutine Algorithm 2 of [7] and the times of the iteration 
that we operate on the subroutine. By simply analyzing the Algorithm, 
we can know that the times is no more than max^if- + -), which can be 
denoted as 0{^). On the other hand, the complexity of the subroutine is 
O(^^m^logn) — km? log n) , see [3] . Hence the complexity of Algorithm 
is 0{{^kw?\ogn). 

By further analyzing the Algorithm, we can easily know that y* is a 

k 

flow on Vi{y*) < h,i = 1,2,- ■ ■ ,k, V{y) < + h*]r],V{y*) < 

i=l 

k 

[{I* - bi) + h*]r], and min < min ^ < l*r]. Hence we only specify 

i=l 

[{I* - 1)(E bi) + {h* - -2r]< V{y*) and min ^ > {I* - l)v - 

i=l 

In terms of the Algorithm, we have e < {I* — l)r] < 1, {h* — l)r] < 

i=l 
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k k 

E and > T^[(E(^* - ^)vbi) + {h* - l)r]]. Therefore 



l+e 

1=1 i=\ 

V{y*) > T^jhi* - + {h* - m 

i=l 

= [{1* - i)(E bi) + {h* - - - k) 
i=i i=i 

+{h* - l)rj] (1) 

> [{1* - b.) + {h* - - 6[(E k) + (E k)] 

1=1 i=l i=l 

> [{l*-l){j:bi) + {h*-l)]v-2v. 



i=l 



Let y/ = E y'iiP),^ = 1, 2, • • • , /c and = E ?/f (^)- Then V^iy*) = 

y/ + since yi^y'i + yf. As 6^ = {I* - l)r]bi, we can have < {I* - l)r]bi 
Suppose Vj{y*) < (/* — l)rjbj — 2rj for some j. Then 



(2) 



< E K + E K° + - 1)^^'. - 2r/. 
On the other hand, E K° < E K° < - 1)^- Hence, (2) imphes 

V{y*) < [{Eil* - i)vbj) + (h* - l)r^] + (I* - l)rjbj - 2n 

(3) 

= [(E(/*-lM) + (/^*-l)^]-2r7. 

It is clearly that (3) is in contradiction with (1). So, Vi{y*) > {I* — 'i-)r]bi — 2r] 
for any i. This implies 



mm 



>(r-l)r/- 



bi min 6j 

This comples the proof. □ 
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Remcirk 2. Recently, Busing and Stiller [1] considered a kind of network flow 
problem arising in line planning. It may be an interesting topic to explore 
the application of the present work in line planning. Moreover, Soleimani- 
damaneh [1 1] investigated the determination of maximal flow in a fuzzy dy- 
namic network with multiple sinks. Mehri [8] studied the inverse maximum 
dynamic flow problem. We believe to extend the present work in a fuzzy 
dynamic network and to address the inverse problem EMCFPSC may also 
be interesting topics. 

6. Concluding remarks 

In this paper, the problem Extended Maximum Concurrent Flow Problem 
with Saturated Capacity is introduced. The existence of its solutions is 
proved. A approximation algorithm to solve the problem is designed, and the 
effectiveness of the algorithm, including the complexity and approximation 
measures, is discussed. The approach we design the algorithm is the speciflc 
contribution of the present work, whose main characters are to construct 
auxiliary network and to implement iterating search through subroutine a 
known algorithm. To solve other network problems with this approach is an 
interesting topic for further researches in the future. 
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for their valuable comments which lead to the improvement of this paper. 
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